Analytic solution for the critical state in superconducting elliptic films 
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A thin superconductor platelet with elliptic shape in a perpendicular magnetic field is considered. 
Using a method originally applied to circular disks, we obtain an approximate analytic solution for 
the two-dimensional critical state of this ellipse. In the limits of the circular disk and the long 
strip this solution is exact, i.e. the current density is constant in the region penetrated by flux. For 
ellipses with arbitrary axis ratio the obtained current density is constant to typically 10~^, and the 
magnetic moment deviates by less than 10~^ from the exact value. This analytic solution is thus 
very accurate. In increasing applied magnetic field, the penetrating flux fronts are approximately 
concentric ellipses whose axis ratio h/a < 1 decreases and shrinks to zero when the flux front reaches 
the center, the long axis staying flnite in the fully penetrated state. Analytic expressions for these 
axes, the sheet current, the magnetic moment, and the perpendicular magnetic fleld are presented 
and discussed. This solution applies also to superconductors with anisotropic critical current if the 
anisotropy has a particular, rather realistic form. 



PACS numbers: 74.60.-w, 74.60.Ge 



I. INTRODUCTION 

Most experiments with high- Tc superconductors deal 
with thin flat samples in a perpendicular magnetic field 
H, for example c-axis-oriented monocristalline platelets 
or films. In this connection, in recent years the problem 
of the critical state of such samples has attracted con- 
siderable interest, see e.g. Ref. || and the references cited 
therein. Exact analytic solutions were obtainecLonly for a 
circular diska and a thin infinitely long strip.ETEl In deriv- 
ing these solutions it was essential that the critical states 
of the disk and the strip have a known symmetry. For 
the strip the critical state distributions of the magnetic 
field and the current do not depend on the longitudinal 
coordinate while the disk has a rotational axis. In both 
cases the directions of circulating currents are fixed and 
known in advance, and thus the critical state equations 
are one-dimensional (ID). In the present paper we obtain 
an approximate analytic solution of a two-dimensional 
(2D) problem, namely, a thin superconductor platelet of 
elliptic shape in a perpendicular magnetic field. This so- 
lution enables one to understand the critical states of real 
superconductors when the above-mentioned symmetry is 
absent. The accuracy of our solution is very high, espe- 
cially for the integral quantities like the magnetization of 
the sample, and it may be used to analyze experimen- 
tal data. Besides this, the analysis of our new solution 
reveals novel qualitative features of the critical state of 
thin superconductors like flux fronts with varying curva- 
ture and the possibility of rotating flux-line arrangments. 

We derive this solution in Sec. II and analyze it in Sec. 
Ill, reproducing the known solutions for disks and strips 
and discussing the interesting new features. In Sec. IV we 



apply our solution to anisotropic superconductors, and in 
Sec. V we summarize the results. 



II. DERIVATION 

Let us place the origin of the coordinate system in the 
center of the sample and let its plane coincide with the 
xy plane. The external magnetic field H is directed along 
the z axis, i.e., along the thickness d of the sample. The 
boundary of the superconductor in the xy plane, F, is 
described by the ellipse (see Fig. 1) 
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where ao and 6o are the semi-axes of the ellipse (ao > 
bf)^ d). The critical current density jc is assumed to be 
constant. The critical-state equations for the thickness- 
integrated current density 3{x,y) = fj^,^j{x,y,z)dz in 
the partly penetrated critical state have the following 
form. In the region between F and the penetrating flux 
front 7 one has 



|J(2^,y)l ^ jcd^ Jc 



(1) 



Here the common assumption has been made that the 
direction of the current does not depend on the coordi- 
nate z across the small thickness of the sample. In the 
flux-free region, i.e. inside 7, the perpendicular induction 
should vanish, 



B,{x,y) = 0. 



(2) 



In addition, the equation 



V-J(x,j/)=0 



(3) 



must hold for any x and y inside the sample. In deriv- 
ing the sokition we follow the method of Mikheenko and 
Kuzovlev.0 When 7 coincides with F, i.e., the flux does 
not penetrate the superconductor, the solution of Eqs. 
(2) and (3) may be found by considering the field on the 
surface of an ideally screening oblate ellipsoid in a uni- 
form magnetic field H\\z. The field at the sm-face of the 
ellipsoid is tangential to it and has the formO 



Ht(a;,2/,z) 



H 



l-N, 



n(zn) 



Here z is the unit vector along z, n is the unit vector nor- 
mal to the surface, and Nzz is the appropriate demagne- 
tizing factor. For a thig oblate ellipsoid with semi-axes 
Qo ^ ^0 3> Co one haaZl 1 — N^z = {co/bQ)E{k) where 
E{k) is the complete elliptic integral of the second kind 
with k^ — 1 ~ h'^/oQ- For example, when aq = 5o (disk) 
one has E{Q) = 7r/2 and 1 — Nzz ~ 7rco/2ao, and when 
ao ^ ^0 (strip) one has E{1) = 1 and 1 — Nzz = co/bo. 
Letting co — *■ and taking into account the relation be- 
tween the surface screening current and the tangential 
field, Jsurf = n X Ht,El we arrive at the sheet current 

J(a;,y) = 2z x Ht , 

where Ht is the field on the upper side of the platelet. 
This result applies to films of any shape; it means that the 
magnitude of the field change across the film thickness 
equals J. In thin ellipses, in the region /(cc/ao, y/bo) < 1 
with f{u, v) = u^ + v'^, one obtains for J = {Jx, Jy). 
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(4) 



where S = 1/^1 - x'^/a^ - yV^o and P ^1- {bo/aof. 
If f{x/ao,y/bo) > 1, we define Fx = Fy = 0. When a 
partially penetrated critical state occurs in the sample, 
we look for the solution of Eqs. (1), (2), (3) as a linear 
combination of distributions like that described by for- 
mulae (4). In other words, we assume that the front of 
flux penetration has the shape of an ellipse with the ratio 
of the semi-axes depending on the depth of the penetra- 
tion, see Fig. 1. This assumption is exact for disks and 
strips and also in the limiting case of a small depth of 
penetration. As will be seen from the results obtained 
below, the deviation from this assumption increases with 
increasing depth but remains very small. The proposed 
approximation thus is sufficient for most practical pur- 
poses. 

Let the penetrating flux front be the ellipse 
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(5) 



with the long axis a depending on the short axis 6. Then, 
according to the above assumption, we write for the com- 
ponents Jx and Jy of the sheet current, 

^.. = -^/V(5')F.„(-^,|;,-^).6', (6) 

where W{b') is a weight function and Fx and Fy are de- 
fined by Eqs. (4). The sheet current (6) with any weight 
W{b') satisfies Eq. (3) due to the definition of the func- 
tions Fx and Fy. In the elliptic region 

this current creates a constant field equal to 



- H / W{b') db' , 

Jb 

which has to compensate the applied field H . Thus, con- 
dition (2) reduces to the equation 



fco 



W{b') db' = 1 



(7) 



The functions W{b) and a (6) must be obtained from the 
equality 



•Jx ' '^y ~ "J c 



(8) 



that follows from Eq. (1). Let x — Q. Then Eq. (8) takes 
the form which was considered in Ref. 0, and we find 



W{b) 



Jc E{k) 



^Hb^i^^wm 



(9) 



with k^ = 1 — [6/a(6)]^. Now we can obtain the function 
Jx{x = 0,y) [while Jy{x = 0, j/) = 0] and calculate the 
total current / circulating in the sample. 



/= / Jx(0, j/)(i2/ = -Jcfcoarccos — 
Jo ^ bo 



(10) 



The H dependence of the penetration depth b{H) along 
the y axis is determined by Eq. (7). This yields. 



'» E{k) db' 



H 



He Jb b'^l-b'^/bl 



(11) 



with He = Jch = Jcrf/TT and fc^ ^ 1 _ [67a(&')]2, see 
Fig. 2. To determine the dependence a(6) let us repeat 
the above reasoning, Eqs. (5) to (10), but now we shall 
consider 6 as a function of a and integrate over a' from a 
to ao- At y = we can find M^(a') and calculate / again, 
obtaining 



f"'" 2 a 

I = — Jy{x, 0) dx — —JcttQ arccos — . 
Jo ■^ Oo 



(12) 



Since the total current crossing any radius must be the 
same, we may equate the currents (10) and (12). This 
yields the dependence of a on fe (Fig. 3), 



a{o) = oo cos — arccos -— 



(13) 



Formulae (4)-(6), (9), (11), and (13) give the solution 
for the critical state in the ellipse. If this solution were 
exact, Eq. (8) would hold for any x and y in the re- 
gion between the curves F and 7. The numerical analysis 
shows that, in reality, the quantity J = ( Jj + JyY^^/Jc 
gradually deviates from 1 as the point {x, y) moves away 
from F. However, this deviation remains small, typically 
10"'' • • • 10~^, and is not visible in Figs. 4-7. In the 
fully or almost fully penetrated state near the points 
X — ±a(& = 0), 2/ = 0, this deviation may reach a 
few percent. Our expressions provide thus a very good 
description of the critical state in elliptic platelets. In 
particular, since this small deviation of J{x, y) from Jc 
occurs mainly near the axis y = 0, its influence on the 
magnetic moment will be negligible, see Sec. III. 

Interestingly, an explicit expression can be obtained 
for the 6-dependence of the magnetic moment m of our 
elliptic disk in the critical state. In the ideal screening 
state any ellipsoid with volume V and semi-axes ag, 60: co 
in a field H\\z has the magnetic moment 



»7iidcai(ao,feo,co) 
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For thin oblate ellipjsoids with oq > 60 3> cq the de- 
magnetizing factor iail Nzz = 1 — {co/bo)E{k) with k^ = 
1 — t'o/'^o- Taking the limit co ^- one obtains the mag- 
netic moment of an ideally screening thin ellipse, 

midcai(ao,M = -f|j|i?. (14) 

For the thin ellipse in the critical state we write 

rbo 

Wib')midcAa{b'),b']db'. (15) 

With the weight (9) the elliptic integral E{k) drops out 
and this integral yields 



2 2 

'-:;JcaobQ 



cos(arcsin6 + e arccos 6) 



cos(arcsin6 — e arccos 6) 

TT~e 



(16) 



where e = bo/ao and b = b/b^. The dependence m{H) is 
obtained from Eqs. (11) and (16), see Fig. 8. 



Next we give the expression for the magnetic field 
Hz(x,y) in the penetrated region between F and 7. Let 
the point (x, y) be located on the ellipse 



X 
"•xy 



y 

xy 



= 1 



(17) 



with semi-axes axy and b^y, where 6 < b^y < &o E^nd 
o-xy — dibxy) satisfies Eq. (13), i.e., the ellipse is a previ- 
ous flux front. Fig. 1. Then we obtain 



Hzix,y)^H, 



, Vl-b'Vbl 



b' K + a^b') E{^,k) 



i{b') V m + b'') 



(18) 



Here £'(1^9, k) is the elliptic integral of the second kind, 

simp^ K/(e + &'^)]^/^ fc2 = 1 - [57a(&')]^ and 



2£,^x^+y^-a 
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l'2-6'2f 



-4(&'2a'2-a;2fo'2-yV' 



1/2 



with a' — a{b') from Eq. (13). The magnetic field out- 
side the superconductor in the z = plane is given by 
the same expression (18) but with the upper boundary 
bxy replaced by 60, and inside the flux front 7 one has 
Hz = 0, see Appendix A. Contour lines and profiles of 
Hz{x,y) are depicted in Figs. 9 - 11. At large distances 
r = (x^ -I- y^)^/^ > flo > &o the field (18) becomes 

iI,(x,y)^i/-TO/(47rr3), 

as expected,Q with magnetic moment m from Eq. (16). 

To end this section we give an elegant representation 
of the sheet current J(a;, y). Substituting Eqs. (4) and 
(9) into Eq. (6), the expression for the components Jx 
and Jy can be written as follows, 



J -J^ 

dy 



Jy 



T ^ 
OX 



gix,y) 



db' 



,^^^T^WM 



'i- 



a'^ib') 



5/2 



(19) 



As in Eq. (18) the boundary bxy is the semi-axis of a pre- 
vious flux front defined by Eq. (17) if the point (x^y) is 
located between the curves F and 7, and bxy = b when 
the point is inside 7 (in the field- free core). For points 
outside the superconductor we put bxy — bo, yielding 
g — 0. The function g{x,y) is the local magnetization 
(or density of current loops) introduced in Ref. Ig; the 
contour lines g{x,y) = const coincide with the current 
stream lines; one has g{x, y) = on the edge F of the 
elliptic sample (and outside), and the volume under the 
"mountain" g(x, y), coincides with the magnetic moment 
(16). We emphasize that the current stream lines in the 
general ellipse are not parallel to the flux front 7, see Fig. 
6, in contrast to the situation in thin disks and strips, and 
in longitudinal geometry. 



III. ANALYSIS 

One easily verifies that the obtained expressions for the 
penetration depth, Eq. (11), the magnetic field, Eq. (18), 
and the currents, Eq|,,(19), go over into the appropriate 
formulae for the diskH when oq = bo and for the stripETEl 
when gq ^ b^. For the strip one has E{k = 1) = 1 and 
formula (11) yields the flux- front position 



(20) 



cosh{H/Hc 



which coincides with Eq. (3.5) of Ref. |j since He = Jc/tt. 
For the disk, E{k = 0) = tt/2 and we arrive at the same 
expression (20) but with He substituted by ttHc/2 — 
Jcf^, see Eq. (19) of Ref. 0. For general ellipses, b{H) 
falls between the curves for the disk and the strip, see Fig. 
2. It should be noted that the penetration along the x 
axis does not go all to the center but a{b) (13) approaches 
the finite value a(0) = aoCOs(7re/2) when 6^0. 

As for the magnetic moment, when ao = 6b, Eq. (16) 
goes over to the appropriate formula for disks | 



2 o/ 1 sinh|/i| 
rridisk = -7- JcOq arccos — - H -2— 

3 V cosh h cosh h 



(21) 



with h = 2H/Jc. In the limit oq ^ bo, Eq. (16) yields 
the magnetic moment of strips with width 2w and length 
2ao, 



mstrip = -2 JcOqw^ tanh{H/Hc) 



(22) 



with an effective width 2w = (2/3) ^/^2&o- This JTistrip 
deviates from the result for a strip with width 2bo by a 
factor 2/3. The reduced effective width is explained as 
follows. A long narrow ellipse with semi-axes ao S> bo 
at each position x may be approximated by a piece of a 
long strip with half width y{x), having a magnetic mo- 
ment proportional to y'^{x)/bo = 1 — x? j a^, cf. Eq. (22). 
Averaging this prefactor over the length 2ao of the ellipse 
one arrives at the reduction factor 2/3. 

The magnetization curves m{H) for disks (21), strips 
(22), and the general eUipse, Eqs. (16) and (11), coin- 
cide almost exactly when they are normalized such that 
their initial slope and saturation value are unity. The 
deviation of these reduced curves Mr{Hr) from the strip 
result tanh(_ffr) are always smaller than 0.013 as shown 
in Fig. 8. The maximum deviation occurs for ellipses 
with excentricity e ~ 60 /ao ~ 0-7. Taking the initial 
slope to'(O) = dm/dH\H=o from Eq. (14), 



m'(0) 



?71idoal 

H 



47r aobl 
' 3 E{k) 



(23) 



(fc^ = 1 — e^, e — bo/cLo), and the saturation value 
fn^at = 'rn{H -^ cxd) from Eq. (16), 



n, n\ 4 2 cos(e7r/2) 

Wsat = mlb ^ 0) = -- JcfloOo — j 2~ 



(24) 



and approximating the normalized curves by Mr{Hr) 
tanh(i?r) with Air = 'm/irisat and Hr — H/ Hi, 



msat J E{k) cos(e7r/2) 



m'(0) 



1-e^ 



(25) 



we get an excellent approximation for the magnetic mo- 
ment of thin ellipses in the Bean model. 



m{H) « TOsat ta,nh{H/Hi) , 



(26) 



with rrisat and Hi from Eqs. (24) and (25). Formula 
(26) is good also for rectangular films, where rrisat is well 
known and m'(0) has to be obtained numerically using 
the method of Ref. S. 

The sheet current J, Eqs. (6Xand (19), reduces fPrfhe 
known results for circular disksH and infinite strips.uQ In 
both limits one obtains the same functional form for the 
magnitudes J{r) = Jcf{r) and J{y) = Jcf{y), with the 
same function: /(y) = 1 for 6 < |y| < 5o and 



I{y) = ~ arctan 

TT 



v/^^^F 



(27) 



for \y\ < b, where c = (1 - &V^§)^^^ = tanh(/i) with 
h = 2H/Jc for the disk and h = ttH/ J^ for the strip. For 
long ellipses with finite length 2ao in these expressions 
bo should be replaced by the effective half width y{x) as 
discussed after Eq. (22). 

There is another known result with which our solution 
may be compared. Namely, in the flux-penetrated re- 
gion where J = Jc is constant, the current stream linss. 
are equidistant and are thus the "distance function" .E2l 
This means all points on a given stream line have the 
same distance from the specimen edge. This property 
applies to thin films of any shape. It applies also in the 
longitudinal geometry, i.e., in both limits of small and 
large thickness. (In the non-penetrated region the cur- 
rent stream lines are different in these two limits; in the 
longitudinal limit one has there j = and J = 0, but 
in the transverse limit J ^ 0). The current stream lines 
in the penetrated region are thus the envelope lines of 
circles with constant radius r [Q < r < bo — b) and with 
center on the specimen edge. Parametrizing our elliptic 
edge by x{ip) — ao cos (p, y{(p) = 60 sin ip, we find for these 
envelopes 



x{ip, r) = (ao - rbo/W) cos ip , 
y{ip, r) = {bo - rao/T4^) sin (p , 



(28) 



with W — {bocos^ ip + aosiii^ ipy^^. The exact stream 
lines (28) have a bend (discontinuity in their slope) in the 
penetrated region on the x axis at a(6) < \x\ < ao — b^/ao, 
but they are smooth at ao — 60/^0 — 1^1 — "^o- Our 
approximate solution coincides with these exact stream 
lines within line thickness. 

From the stream lines in the critical state, Eq. (28), 
one obtains the exact magnetic moment of thin ellipses 
in the fully penetrated critical state, 



<aT* = -Pcaobl [E{k) - '-K{k) ] , (29) 

where K{k) and E{k) are the elhptic integrals of the 
first and second kind with k^ = f — e^, e — ba/ao, see 
Appendix A. Comparing the exact result (29) with our 
approximate rrigatj Eq. (24), we find that the deviation is 
extremely small: 

Writing /i(e) = cos(e7r/2)/(f - e^) and /2(e) = 
E{k) - {e^/2)K{k) we find /i(f) = /2(f) = 7r/4 (disk), 
/i(0) = /2(0) = f (strip), and < /i - /2 < 0.001317. 
The maximum deviation /i — /2 — 0.001317 occurs at 
e — 0.5043. The deviation is an almost symmetric func- 
tion of the excentricity e, /i(e)— /2(e) ~ 0.00138 sin^(e7r). 
With Eq. (29) the field Hi (25) which enters formula (26) 
takes the exact value 



HI' 



to'(O) 



H,E{k)[Eik)~-K{k)]. (30) 



The flux fronts, current stream lines, and contour lines 
of the magnetic field Hz{x,y) in the general elliptic film 
do not coincide, as opposed to the situation in circular 
disks, long strips, and in longitudinal geometry. These 
features of the critical state in non-circular thin super- 
conductors of finite size arej-.Sjeen even more clearly in 
square and rectangular filmsJjiHI but while those geome- 
tries were solved numerically, we now have an analytic 
solution for the elliptic film . 

Another novel feature of the critical state in thin el- 
liptic films as compared to the critical state in circu- 
lar disks and infinite strips, is that the direction of the 
sheet current J at a given point {x,y) changes during 
the penetration of flux, see Fig. 7. This fact may have 
an important physical consequence for films or platelets 
with thickness exceeding the London penetration depth 
A. As is known, tdTli in flat superconductors with d:$> X 
in moderate fields H there is a flux and current-free 
core described by \z\ < Zcoro(a;,y). In the cases of 
strips and disks this core was depicted, e.g., in Refs. 
|l^,|l6|. Critical currents circulate only outside the core, 
for Zcorc(x,y) < l^l < d/2. In this shell the critical state 
has the usual Bean form, but the field gradient occurs 
across the thickness of the sample since this state is forced 
by the screening currents, and the flux lines are almost 
parallel to the flat surface. The outer rim (equator) of the 
flux-free core coincides with the penetrating flux front 
since in the region inside the flux front the perpendicular 
magnetic field Bz is practically zero. 

Now, since the tangential field at the flat surface is per- 
pendicular to the sheet current J, the changing direction 
of 3(x,y;b) during flux penetration (i.e. with increasing 
H and decreasing b) means that the direction of the pen- 
etrating flux lines in the region below 7 also changes. 
Thus, flux lines with gradually rotating orientation enter 
through the flat surface at each point away from the sym- 
metry axes. These U-shaped flux lines move towards the 
specimen center, dragged by the Lorentz force, which is 
balanced by bulk pinning. As a result, in the region inside 



7 the layers of the flux- line lattice exhibit torsion relative 
to each other and therefore carry a longitudinal current 
component, i.e., the local current density is not exactly 
perpendicular to the flux lines. In principlft,4;bis may lead 
to instabilities and flux-cutting processes.tZnlj That sit- 
uation may be related to the penetration of a rotatina 
field component discussed by BeanE2l and by Gilchrist .Ell 
In deriving the solution we have assunied|-t]hat B — 
/ioH and the so called geometrical barrierc3~E3 is negli- 
gible. These assumptions hold if the characteristic mag- 
netic field jcd considerably exceeds the lower critical field 
Hci ■ In the case jcd < Hd and H < Hd the criticaLstate 
problem becpHies more complicated. AnalyticalE3nJ and 
numericalcj'Ej methods to analyze this general problem 
were recently proposed. In particular, if our assumption 
B = fioH is replaced by the coa£ct H{B) (obtained, e.g., 
from Ginzburg-Landau theorytj) then flux does not pen- 
etrate until the applied field H has reached the field of 
first penetration of flux. Hen ~ Hd tanh ^Jcd/2w, where 
c « 0.36 for strips mith half width w and c « 0.67 for 
disks with radius w\ 



IV. ANISOTROPIC CRITICAL CURRENT 

The obtained solution allows one to analyze the crit- 
ical state of fiat superconductors when the critical cur- 
rent density jc is anisotropic in the plane of the film. 
Such anisotropy can occur, for example, if there are twin 
boundaries or other extended or inclined defects, or if 
the superconductor itself is anisotropic. A controlled in- 
plane anisotropy of jc of more complicated nature may 
be induced by applying a magnetic field parallel to the 
filmEj To explain our analysis, we perform a transforma- 
tion of coordinates: 



ax', y = (3y', 



(31) 



where a and /3 are some constants and the prime denotes 
the transformed quantities. Under this transformation 
the ellipse with semi-axes oq and bo goes over to an el- 
lipse with semi-axes Og — ao/ct, 6q — bo/P, while any 
small line element, dl, directed at an angle (p relative to 
the X axis, is transformed into the element 



dl' = dl 



? -2x1/2 

COS"" ip sm if ^ 



/?2 



which is directed at an angle if' with 

, a 

tan (f = — tan ip . 

P 

Since the current crossing the line element must be in- 
variant, we obtain the following transformation of the 
sheet current flowing at an angle ijj — ip± 7r/2 relative to 
the X axis: 



J'(V'') = JW(3 



sm ip 



1/2 



Thus, if the anisotropic critical current density in some 
elhptic sample with semi- axes Aq, 6q can be approximated 
by the expression 



J'^i^') =const•(l + (5sin2V'')^^^ 



(32) 



in which 6 is some constant (S > —1), we may reduce the 
critical state problem of this anisotropic sample to the 
isotropic problem by performing the transformation (31) 
from x' , y' to x, y with (3—1 and a — \/l + 5. After 
the transformation we shall have an elliptic sample with 
semi-axes ag — vT+Jag, &o = ^o ^^-^ with isotropic 
critical current Jc — const, and we can use the solution 
obtained in Sec. II. 

Thus, our solution permits to analyze (at least quali- 
tatively) the features of the critical state of anisotropic 
superconducting films. For example, for a circular disk 
with radius R and with anisotropic critical current obey- 
ing Eq. (32) with (5 = 1, one obtains the transformed 
ellipse with semi-axes ag = v^i?, bo = R and the cur- 
rent distribution and stream lines depicted in Figs. 4 - 
7 for e = 0.7, and the magnetic field depicted in Figs. 
9 - 10. These pictures then are transformed back to 
the coordinates of the original circular disk by shrink- 
ing them along a; by a factor of 0.7. Figure 12 shows 
another example of a circular disk with larger anisotropy 
^ = 3, yicldinff Jc(7r/2)/Jc(0) = ao/6o = 2. The re- 
sults of Ref. g^ for square and rectangular plates with 
anisotropic critical current are in agreement with this ap- 
proach. Our solution may also be used to include thfi-ip. 
plane anisotropy of the critical current into theorieaU'ES 
which analyze the critical state in flat superconductors 
taking into account the current anisotropy (B/B depen- 
dence) associated with the curvature of the flux lines. 



V. CONCLUSION 

In this paper we derived an approximate analytic so- 
lution of the Bean critical state model for thin elliptic 
superconductors with semi-axes ag and bo = eag in a 
perpendicular field H. The accuracy of this solution 
is high and the known limiting cases of a circular disk 
and a long strip are recovered exactly. With increasing 
H the penetrating flux front is an ellipse with axes a 
and b related by Eq. (13) such that at full penetration, 
when 6 = 0, a{b) stays finite, i.e., there is a section of 
length 2a(0) = 2ao cos(7re/2) on the x axis where regions 
with different orientation of the curreiijLjruaet, similar to 
the situation in rectangular platelets. ElO'EZl The current 
stream lines in the penetrated region deviate from ellipses 
and, in the fully penetrated state, they have a sharp bend 
on the X axis in the interval |x| < gq — b'^/ao that is 
slightly wider than 2a(0). 

Our analytic solution shows that in non-circular finite 
superconductor films at any point not located near a sym- 
metry axis or near the edge, when H is increased the di- 
rection of the circulating sheet current changes until the 



flux front passes through this point, see Fig. 7. There- 
fore, in samples with thickness d exceeding the London 
penetration depth A, the orientation of the flux lines pen- 
etrating from the two flat surfaces changes also, since at 
the flat surfaces the flux lines are at a right angle to the 
sheet current. As a consequence, the flux-line arrang- 
ment rotates along z in the region where B^ = 0. Fur- 
thermore, it cannot be ruled out that also in the region 
with Bz ^ Q the flux lines are twisted, i.e., their cur- 
vature is not planar. If this twist occurs, it means that 
the orientation of the current at flxed x, y but different 
depth z varies, and the integrated current density J{x, y) 
is not necessarily equal to jcd or exactly constant even 
in the flux-penetrated region. The precise calculation 
of the flux and current distributions inside non-circular 
disks even of small thickness d > X then becomes an in- 
tricate 3D problem. The possibility of twisted flux lines 
in the penetrated region of non-circular thin disks should 
thus be considered in more detail in future work. 
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APPENDIX A: EVALUATION 

The integrals (11) for H/H^, (18) for Hz{x,y), and 
(19) for Jx, Jy, and g{x,y), like the elliptic integral of 
the second kind, 



E{ip,k) ^ E{s,m) 



ms 



2^2 



1 



2^2 



dt 



(Al) 



with s = sin ip, m = k^ , may be evaluated with high pre- 
cision by a substitution of variables. We flrst transform 
the integrals over b' into integrals over a variable t in the 
interval < t < 1, like Eq. (Al). Noting that these in- 
tegrands may have poles ex 1/Vi and ex l/Vl — t at the 
boundaries, we use a substitution 



1 .1 

f{t)dt^ / f[t{u)]t'{u)du (A2) 



which has a weight function t'(u) = dt/du that vanishes 
with a high power of u and 1 — uaXu — t = Q and at 
M = t = 1. A good such choice is 

t = g{u) = 'ibu^ - Mu^ + 7Qu^ - 20u^ , 

t' = g'{u) = 140u^(l - uf . (A3) 

An even better choice is to iterate (A3) once, writing 

t^g[g{u)], t'^g'{u)g'[g{u)]. (A4) 

With (A4) one has t ex u^^ and t' oc u^^ near t = 0, and 
1 - t ex (1 - u)i6 and t' ex (1 - u)'^^ near t = 1. The pole 



f{t) oc 1/Vt in the original integrand is now removed in 
the substituted integrand F{u) = f[t{u)]t'{u) ex u^ and 
similarly near t = I. Since the new integrand vanishes 
rapidly at the boundaries of the integral, one may use an 
integration method with constant weights, 



.1 , N 

I^ F{u)du^-J2P 



i- 1/2 

N 



(A5) 



This method achieves high precision « 10^^ even with a 
small grid number A^ = 20 ... 30. 

In Eqs. (11), (14), etc., E{k) is given by the definition 
(Al) with s — 1, and for K{k) in Eq. (29) a similar defi- 
nition applies. All these elliptic integrals are easily com- 
puted with high accuracy by this integration method. 

Before the integrals for Hz{x,y) and g{x,y) can be 
evaluated for a 2D grid of points {x,y), one has to find 
the elliptic flux front that passed through the point {x,y). 
This may be done by finding the zeros of the function 



Pibxy 



X 

'a?~ 

"'xy 



V 



1 



with axy = cos[6oarccos(&a;iy/6o)], cf. Eqs. (13) and (17) 
with length unit ao = 1. This is achieved by starting 
with bxy = b (the short semi-axis of the inner flux front) 
and iterating 6a;j, <— bxy-2ep{hxy)/[p{bxy + e)-p{hxy-e)] 
(e <C 1) a few times. A similar procedure may be used 
to invert the relation H{b), Eq. (11), to obtain b{H). 
To obtain the correct g{x,y) also inside the flux front, 
and Hz{x,y) also outside the superconductor, we chose 
bxy — bo for all points outside the ellipse with semi-axes 
ao, bo, and bxy = b for all points inside the ellipse with 
semi-axes a, b by putting bxy ^- max[&, rmn{bxy, bo)]- 
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FIG. 1. The three ellipses defining the edge F of the film 
with semi-axes ao, feo; an earlier fiux front passing through a 
given point {x,y) (semi-axes ttxy, bxy, dashed line); and the 
present flux front 7 with semi-axes a, b. 
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FIG. 2. The semi-axes a and h of the elhptic flux fronts 
plotted as reduced penetration depths (ao — a) /bo (dashed 
hues) and {bo — b)/bo (solid lines) versus the applied field H 
(in units of He = Jc/tt) for thin ellipses with excentricities 
e = bo/ao = 1, 0.9, 0.7, 0.5, 0.4, 0.3, 0.2, 0.1, 0.05, 0. The 
curve for the disk [e = 1, a — b, Eq. (20) with He replaced by 
Jc/2] is marked by circles. For strips, b [Eq. (20)] is marked 
by crosses, and one has ao — a = 0. 
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FIG. 4. 3D plot of the magnitude of the sheet current 
\3{x,y)\, Eq. (19), in a thin ellipse with bo/ao = 0.7 at 
b/ao = 0.2. In the flux-penetrated region the plotted J/Jc 
deviates from unity by less than 10~^. 
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FIG. 5. Profiles of the sheet current Jx{0,y) in a thin el- 
lipse with bo/ao = 0.7 for various degrees of flux penetration 
b/ao = 0.69, 0.6, 0.5, . . . , 0.1, 0.05 (from right to left). Note 
the precise saturation to J = Jc, which demonstrates the ac- 
curacy of our approximation. 



FIG. 3. The long semi-axis a as function of the short 
semi-axis b of the penetrating elliptic flux front for various 
excentricities e = bo/ao — 1, 0.9, 0.8, . . . , 0.1, 0, Eq. (13). 






FIG. 6. Stream lines (solid lines) of the sheet current 
i}{x,y) in a thin ellipse with (a) bo/ao — 0.7, b/ao = 0.2 
(yielding a = 0.624) and (b) bo/ao = 0.5, b/ao =0.1 (yielding 
a — 0.775). Also shown are five previous flux fronts (dashed 
lines) at equidistant values of b. The edge of the film and 
the inner flux front are indicated by circles and a bold line. 
Note that in general the sheet current does not flow parallel 
to the flux fronts and that its magnitude is constant in the 
flux-penetrated region. 



FIG. 7. The stream lines of the sheet current in the same 
ellipse {bo/ao = 0.7) at two different times during flux pen- 
etration: at b/ao ~ 0.65 (solid lines, almost ideal screen- 
ing) and b/ao ~ 0.01 (dashed lines, almost full penetration). 
Note that the orientation of the sheet current at some places 
changes strongly during the penetration, leading to a rotating 
vortex arrangmont across the thickness of the film. 
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FIG. 8. The magnetic moment (magnetization curve) of 
thin ellipses in a perpendicular magnetic field in the Bean 
model, plotted in normalized form as deviation from the strip 
result Mr = ta.nh{Hr) versus a reduced field Hr = H / H\ for 
excentricities e = bo/ao = 1, 0.9, 0.7, 0.5, 0.4, 0.3, 0.2, 0.1, 
and 0.05, see Eqs. (11), (16), and (24) - (26). 




FIG. 9. Contour lines (solid lines) of the perpendicular 
magnetic field Hz{x,y) of a thin ellipse with bo/ao — 0.7 in 
the plane z = for partial flux penetration with b/ao — 0.2, 
as in Fig. 6a. Note that these contours are nearly parallel to 
the elliptic flux fronts shown as dashed lines for b/ao = 0.7, 
0.6, 0.5, 0.4, 0.3, and 0.2. The edge of the specimen and the 
inner flux front are marked by bold lines. Inside the inner flux 
front one has exactly Hz{x,y) — 0. 




FIG. 10. Profiles of the magnetic field H^{x,y) of Fig. 9 
plotted along y (left) and along x (right) for discrete equidis- 
tant values of x or y, respectively. The dashed line marks the 
apphed field H/Hc = 2.39. 




FIG. 11. 3D plot of the magnetic field Hz{x,y) for 
bo/ao — 0.7 and b/ao = 0.2 as in Figs. 9 and 10. The logarith- 
mic infinity at the edge of the ellipse is cut off at Hz/Hc = 3.5 
to limit the irregular peaks caused by the equidistant grid used 
for this plot. 
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FIG. 12. Current stream lines (solid curves) and flux fronts 
(dashed curves) in a circular disk with anisotropic critical cur- 
rent density, Jcy/Jc^^ = J{tv/2)/Jc{0) = 2 [5 = 3 in Eq. (32)] 
for half {b/bo = 0.5, left) and fuU {b/bo = 0.02, right) pen- 
etration of flux. Shown at the bottom is the isotropic ellipse 
with axis ratio bo/ao = 0.5 to which this problem transforms. 
This ellipse exhibits equidistant current stream lines outside 
the inner flux front. 
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